Complete set of operational measures for the characterization of 3— qubit entanglement 
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We characterize the entanglement contained in a pure three-qubit state via operational entangle- 
ment measures. To this end we derive a new decomposition for arbitrary 3-qubit states which is 
characterized by five parameters (up to local unitary operations). We show that these parameters 
are uniquely determined by bipartite entanglement measures. These quantities measure the entan- 
glement required to generate the state following a particular preparation procedure and have a clear 
physical meaning. Moreover, we show that the classification of states obtained in this way is strongly 
related to the one obtained when considering general local operations and classical communication. 
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Entanglement is at the core of many of the applications 
of quantum information theory and quantum computa- 
tion and plays a key role in the foundations of quantum 
mechanics. Therefore, a great amount of theoretical ef- 
fort has been performed in recent years to grasp this phe- 
nomenon, in particular regarding its characterization and 
quantification as well as its convertibility properties [![. 
Whereas bipartite entanglement is well understood, mul- 
tipartite entanglement is much more subtle. In fact, our 
understanding of the nonlocal properties of many-body 
states is far from complete even in the simplest case of 
just three subsystems. Our knowledge of bipartite entan- 
glement stems from the fact that in the asymptotic limit 
of many copies of any given state, there is a unique op- 
timal rate at which it can be reversibly transformed into 
the maximally entangled state 0. However, such an ap- 
proach seems formidable in the multipartite regime [3|. 
Nevertheless, different classes of multipartite entangled 
states have been identified according to their convert- 
ibility properties [3| and several entanglement measures 
have been proposed like the tangle Q , the localizable en- 
tanglement , and the entanglement of assistance 0, @] 
to cite a few. A fundamental property of entanglement 
is that it is invariant under local unitary (LU) opera- 
tions. This has also led to the study of complete sets 
of (polynomial) invariants under this kind of operations 
[§] and the necessary and sufficient conditions for LU- 
equi valence have recently been provided [ltj. However, 
a complete classification of LU classes with operationally 
meaningful measures was still lacking. In this paper we 
solve this problem for the case of three qubits. That is, 
we provide a complete set of operational entanglement 
measures which characterize uniquely all 3-qubit states 
with the same entanglement properties. These measures, 
which are easily computable, characterize the different 
forms of bipartite entanglement involved in the genera- 
tion of the state following a particular preparation pro- 
cedure. Our results provide a physical classification of 
pure multipartite entanglement and, hopefully, will pave 
the way for new applications of many-body states in the 
light of quantum information theory. 

After removing all the free parameters due to LU oper- 



ations, two-qubit pure states have one non-local param- 
eter: the Schmidt coefficient. It captures all the infor- 
mation about the entanglement of the state. For 3-qubit 
pure states the number of nonlocal parameters is five 
UJ. Thus, one would expect five measures to charac- 
terize entanglement in this case. Surprisingly, one needs 
to consider a set of six LU— invariants to characterize the 
LU-equivalence classes 0. This is because of the coun- 
terintuitive fact that there exist states which are not LU- 
equivalent to its complex conjugate [5, 10] (which will be 
taken with respect to the computational basis and will 
be denoted by *). The sixth invariant is required to dis- 
criminate these states. The class of states (not) being 
LU-equivalent to its complex conjugate will be referred 
to as CLU (NCLU). The existence of NCLU states is a 
striking feature of multipartite entanglement which does 
not exist in the bipartite setting. Since complex conjuga- 
tion corresponds to the redefinition of the complex unit, 
and \ip*) will have the same entanglement properties 
10j . Thus, one cannot expect operationally meaningful 
measures to discriminate between these two states. How- 
ever, as we will show, it is possible to identify states as 
CLU or NCLU by operational tasks. 

The outline of the remainder of the paper is the follow- 
ing. First, we review some known results and introduce 
our notation. Then, we derive a decomposition of 3-qubit 
states, which is characterized by five parameters. This 
form leads us naturally to a physical process generating 
arbitrary 3-qubit states. Next, we characterize this pro- 
cess by five operational entanglement measures and show 
that these five measures together with one binary mea- 
sure characterize the LU-equivalence classes of 3-qubit 
states (up to complex conjugation). Moreover, we show 
that the CLU-class can be identified by the values these 
measures take. Last, we show that the characterization 
of the LU-equivalence classes within CLU is strongly re- 
lated to the characterization obtained by considering lo- 
cal operations and classical communication (LOCC). 

Bipartite entanglement will be measured by the Von 
Neumann entropy of the reduced state, E(\\jj) 12) — 
S(pi), and its convex roof extension for mixed 
states, the entanglement of formation E(p) — 
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min K,l^}} Ej KjEifa) for p = £\ ^j\i>j){^j\- For two- 
qubit states, E = (£(C) is a simple monotonously in- 
creasing function of the so called concurrence, which is 
defined as C(ip) = \(tp\tp}\, where = a y ® cr y \ip*) and 
a x,y,z denote the Pauli operators, and with the convex 
roof construction for mixed states [l[ . An important con- 
cept used here is entanglement of assistance. One of the 
parties (say 1) assists the other two (2 and 3) in obtain- 
ing a particular bipartite entangled state by performing 
measurements on his particle (after possibly adding aux- 
iliary qubits). Let \ip) = ^|0)i |^o)23 + V 1 — _p|l)i|^i>23 
be the Schmidt decomposition of \ip) in the 1 1 23 splitting 
[Tjj . The states, \xi), 23 obtain with probability qt after 
a measurement in 1 are then given by 

1 

3=0 

where U'U = 1. Here, i = 1, . , . , d, where d equals the 
effective dimension of the subspace in which 1 measures 
after the (possible) addition of ancillas. Thus, the states, 
\xi), 2 and 3 can get in this way are given by all ensemble 
decompositions of p 23 , {ftj, !</>,) }, 7 r, yielding the proba- 
bility of obtaining outcome j [141 ] . Hence, the entan- 
glement of assistance, which is defined as the maximum 
entanglement on average one can generate by this proce- 
dure 0, is then given by E% 3 = max {Trj ^. )} £\ TTjE^j). 
Let us now introduce the symmetric matrix t, t^j — 
y /PiPj(4'i\' , Pj}j and the shorthand notation Ci = (ipi\ipi) 
and c = (^o|V'i)- Unless stated otherwise we will choose 
the global phases of \ipo) and \ipi) such that both Co and 
ci are nonnegative. There exist closed formulas for both 
the concurrence C{p 23 ) and the concurrence of assistance 
C a {p 23 ) in terms of the singular values of r [2,13]. 

We use now the notion of assistance to show that any 
3-qubit state can be written as an equal superposition 
of biseparable states which contain the same amount of 
entanglement fljj . I. e. it can be written up to LUs as 

|V) = -^(|0)i|V.)as + |l>iCa ® W.>as). (2) 

Here, \ip s ) = a\QQ) + b\U) with a > b has Schmidt de- 
composition and U 2 = Z(a)Y(J3)Z(-y) and U 3 = Y{fi') 
with Y, Z(£) = e l £ ay >< s respectively. Thus, any state is 
characterized by the five parameters, {E(ip s ), a, /3, 7, /?'}, 
where E{$„) £ [£(C 23 ), £(C 2 a 3 )], a, 7 € \-tc/2,w/2] and 
/3,/3' G [0,tt/2] (see Appendix A). 

To see that any state can be written as in Eq. ([2]) 
(up to LUs), we consider a measurement on qubit 1 
such that both outcomes are equally likely, i.e. q\ = 
q 2 = 1/2 and E{\ Xl )) = E(\x 2 )) (see Eq. ©) 
Since |^o)-L|V'i) we find that the unitary in Eq. (Q} 
can be taken without loss of generality of the form 



|pc + (l-p)cie 2 ^±2 v /p(l - p)ce lul \ and, thus, E(\x )) = 
E(\ X1 }) iff 

ui = arctan [ - ^~ ; CQ ^ ar gg ) (3) 

\pc - (1 - p)ci J 

Since this equation has a solution for any given value of 
argc, this proves the statement. The particular form of 
the unitaries, U 2 ,U 3 , is achieved by taking the Euler ZYZ 
decomposition and using the fact that Z(£) <g) t\ip s ) — 
1® Z(^)\ijj s }. Notice that the decomposition @ (i. e. the 
choice of {E(ip s ), a, f3, 7, /?'}) can be made unique (see 
Appendix A). 

The decomposition ([2]) shows that any 3-qubit state 
can be generated in the following way (see Fig. 1): Qubit 
1 is initialized in the state |+) and qubits 2 and 3 are 
prepared in the state \ip s ) [21] . Then qubit 1 is attached 
to qubit 2 via the controlled nonlocal unitary U^ 2 and to 
qubit 3 via U} 3 where U\ l = |0> 1 <0| ® % + |l)i(l| ® U h 
i. e. |V) = C7 c 13 f/ C 12 |+) 1 |V S )23. Note that [t/ c 12 , C/ c 13 ] = 0. 



a) b) 
3 k> 2 3 k) 2 




Figure 1: a) The decomposition (|12|) of an arbitrary 3-qubit 
state: \ip) = Ul 3 Ul 2 \+)i\ r ips)23 and b) the generation process: 
each party makes a Bell measurement on its qubit and on one 
of its share of the corresponding CJ state 



Let us now briefly discuss a method to generate the 
state using entanglement. It is a well-known fact that 
any nonlocal map, £ , corresponding to a state \^s) via 
the Choi Jamiolkowski (CJ) isomorphism can be (proba- 
bilistically) implemented using a system prepared in the 
state \^e) and local Bell-measurements [3j (see Fig. 1). 
Since the implementation of each control gate, t/^ 2 and 
L/ c 13 , requires four local Bell-measurements, we obtain, 
depending on the measurement result, one out of the 
4 4 different states, {U^a^ <g> af 2 U^ 2 a} 3 <g> cr^ | +) 1 |-0 S ) 23 } 
where <tq — 1. Due to the symmetry of these states, it 
can be shown that this set coincides up to LUs to the 
set <*></, = {Ul 3 Ul 2 o 2 \+) 11^)23}, containing only four 
states (see Appendix B). Each state in this set occurs 
with equal probability. Note that if the aim was to gen- 
erate any other state within S^, one would again, de- 
pending on the measurement outcome, end up with one 
of the states in That is, the set is closed under 
this generation process. In the following we derive bi- 
partite entanglement measures, which characterize such 
a generation process completely. This will finally lead to 
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our main result: the operational characterization of the 
entanglement contained in 3-qubit states. 

In order to characterize this process we need on the 
one hand E(\ip s )) and on the other we need to charac- 
terize the entanglement properties of 2-qubit controlled 
unitaries. Recall that any 2-qubit unitary can be written 
as Vi ® V 2 U n iVs ® Vi [ltjj . Here, the nonlocal part U n i 
is diagonal in the Bell basis and depends in general on 
three parameters. Let us consider the control gate [7 12 
and write without loss of generality XI % — WZ(2 , y)W\ 
where W — Z(a)Y(b) with a, b <E R. Then, in this case 
it can be easily checked that the nonlocal content U n i of 
U^ 2 depends on a single parameter, 7, which is in turn 
uniquely characterized by Ei mp (U^ 2 ) , which denotes the 
amount of entanglement (in the splitting 1|2) contained 
in the CJ-state corresponding to U^ 2 , i. e. the amount of 
entanglement needed to implement U\ 2 ■ Thus, knowing 
just U n i is not enough to characterize J7 C 12 as this does 
not specify the other parameters a and b, which depend 
on the LUs {Vi}. These are as well relevant to character- 
ize the entanglement properties of U^ 2 as can be easily 
understood as follows. Let U}^ 2 = |0)<0|(8)l+|l}(l|(g)f/ 2 
and U}^ 2 = \® |0> <0| + U 2 ® 1 1 > < 1 1 (i. c. the same con- 
trolled unitary but with different control qubit). It can 
be seen that the nonlocal parts of these operations are 
equivalent. Thus, the choice of the local unitaries Vi 
determines which qubit is the control qubit: U}^ 2 = 
X®WZ{~f)U n il®W^ and U^ 2 = WZft) ®U7njW+®l. 
Obviously the information about the control qubit is im- 
portant, e.g. the CNOT gate acting on | + 0), generates 
either a product state or maximally entangled state de- 
pending on which qubit is the control qubit. In summary, 
Ei m p alone does not characterize the entanglement prop- 
erties of a controlled unitary. One needs further measures 
to characterize the local unitaries Vi- This will allow to 
specify U^ 2 uniquely. As we will see, this will make it also 
possible to characterize our generation process and will 
lead to a complete classification of 3-qubit entanglement. 
The most natural choice is then the entanglement created 
after the implementation of the controlled unitary 



Since, in the considered case, the controlling direction is 
not specified by U n i a reasonable choice is the entangle- 



where EJ ain {Ul 2 ) {E^ am (Ul 2 )) measures the amount of 
entanglement generated between system 1 and systems 
2 and 3 by applying the control gate C/ c 12 to the in- 
put state |+)i|'0s)23 and the direction of the arrow in- 
dicates whether system 1 (— >) or system 2 («— ) is the 
control qubit. To see that this is indeed the case, we 
need to verify that {E(\tp s )),a, /3, 7,/?'} in Eq. ([2]) can 
be obtained from the values of {Ei}. Since for all 
the bipartite entanglement measures at least one sub- 
system belongs to an effective two-dimensional space, 
they are in one-to-one correspondence with the concur- 
rence, and therefore with the trace of the squared re- 
duced density matrix. Thus, E3 is characterized by the 
value of (1 + \(<l>+\Y(fi') ® l|0+)| 2 )/2 = (1 + cos 2 p')/2, 
which shows the one-to-one correspondence of E3 and 
0' (recall that 0' G [0,tt/2]). Analogously, E 2 , E A 
and E5 are respectively in unique correspondence with 
cos 2 ^cos 2 (a + 7), cos 2 ^[(a 2 - b 2 ) 2 + 4a 2 & 2 cos 2 (a + 7)] 
and a 4 + b 4 + 2a 2 b 2 [cos 2 (a + 7) cos 2 + sin 2 (a - 7) sin 2 0] . 
Clearly then, E 2 and E4 determine uniquely and a + 7 
modulo 7r, which has only (at most) four LU inequivalent 
solutions which are compatible with £5: {a, 7}, {7, a}, 
{a±7r/2, 7±7r/2} and {7±7r/2, a±7r/2}, where the signs 
in the last two cases are chosen such that the phases are 
in the desired region. This set of states coincides with 
S^p [231 ] . Hence, a one-to-one correspondence is made be- 
tween the entanglement measures {Ei} and the set of 
states obtained through our generation process. 

Let us now have a closer look at the set 
= {f/ c 13 f/ c 12 a 2 |+) 1 |^ s ) 23 }. Note that = 
W), W*)h where \f) = C/ C 13 [/ C 12 a 2 |+) 1 |^) 2 3. 
As discussed in the introduction there is no operational 
measure to distinguish a state from its complex conju- 
gate, hence it is not surprising to have this degener- 
acy. On the contrary, one could ask if there is an op- 
erational procedure to further distinguish \ip) and l^'). 
This extra bit of information can be easily obtained by 
considering the value of the final entanglement gener- 
ated in the splitting 1|23, £'i|23(V')- It can be seen that 
13. E i\m{i>) = #i|23(V>') iff \i>) or \tp*) are LU-equivalent 
to \i>') (see Appendix C). Thus, define the binary mea- 



sure Eq 



E^ as E 6 



(Eg = 1) when £i| 23 takes 



ment gained with the two possibilities, £w n (£^r ) and the minimal (maximal) possible value inside S^. The 



E'gairSVc )i f° r some physically motivated choice of input 
state. It is easy to see that these measures will specify a 
and b. 

Now we are in the position to state our main result, 
the operational characterization of the LU classes up to 
complex conjugation of 3-qubit states. First, we have 
that the set S^p, that is the generation process described 
above, is uniquely characterized by the values of the fol- 
lowing five bipartite operational measures 22 1 



Ei — E(ifj s ), E 2 — E imp (U} 2 ), E3 — E imp (U^ 3 ) 



E 4 



E gtin( U c )> E 5 = E min( U < 



12\ 



gain * 



(4) 



operational measures {Ei}® =1 characterize uniquely (up 
to complex conjugation) the LU classes of 3-qubit states, 
and, therefore, their entanglement properties. 

Despite the fact that Ei(ip) = Ei(tp*) Vi, as expected, 
our approach can, nevertheless, provide also information 
on this issue: the values the measures {Ei} take depend 
on whether \ip) belongs to the CLU or to the NCLU class. 
More explicitly, we have that <E CLU iff either E\ = 
£(C 23 ) or Ei = (£(C 2 a 3 ) (while £(C 23 ) < E x < <£(C 2 Q 3 ) 
for NCLU states) (see Appendix D). In other words, it 
takes the same bipartite entanglement resources to gen- 
erate a state and its complex conjugate, but while for 
CLU states the isolated qubit is attached to some bipar- 
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tite state whose entanglement is either the maximal or 
minimal possible, the NCLU states require some inter- 
mediate amount of entanglement. 

Furthermore, a refined classification of the CLU class 
can be achieved by considering the value of E± . We have 
the following subclasses of states \ip) with E\ — Ei{ip) 
(in what follows c€l and 4>i <E R 2 ) (see Appendix E): 

1. E 1 = £(C 23 ) = <£(C 2 a 3 ) iff |V>) G W-class 0. 

2. Ek = <£(C 2 a 3 ) ± £(C 23 ) iff W (x |000) + clfafafa). 

3. E 1 = <£(C 23 ) £{Ch) iff |V^> oc |000)+e lc |<^ 2( /) 3 ). 

4. 3|^ s ), forEq. dHJ) such that E^}) = <S(C 23 ) 
A EkiW,)) = <£(C 2 ° 3 ) (<£(C7 23 ) ? (£(C&)) iff |V) cx 

(|000> + |010203>)- 

This can be regarded as a classification beyond the 
stochastic LOCC paradigm Q since the last three classes 
correspond to GHZ-type states. This classification is re- 
lated to the classification of states according to determin- 
istic LOCC operations. It follows from the results of [3] 
that classes I, 2 and 3 are closed under these transfor- 
mations, i. e. deterministic LOCC can only take a state 
to another state within the same class. Class 4 is com- 
posed of those states for which \ip s ) in the decomposition 
given by Eq. ((2]) is not unique. In fact, for these states 
Ei can take any value between <£(C 23 ) and <£(C 23 ) (cf. 
Appendix E). This richness is directly connected to the 
LOCC properties of this class of states: these are the only 
ones that might be transformed to other classes [3] (in- 
cluding NCLU states). Note that, due to the symmetry of 
the states in any class, this classification is invariant un- 
der the considered splitting, that is, if e.g. E\ — £(C 23 ) 
in splitting 1 | 23, it is so in any other splitting. 

In conclusion, we derived a decomposition for 3-qubit 
pure states, which, as the Schmidt decomposition for bi- 
partite states, can be easily computed. This decompo- 
sition leads naturally to a generation process, which is 
characterized by the bipartite entanglement measures, 
{Ei}f =1 (Eq. Q). The set {£jf =1 together with the 
binary measure Eq forms a complete set of operational 
measures identifying the different LU-equivalence classes 
(up to complex conjugation). Hence, the nonlocal prop- 
erties of 3-qubits states have been operationally char- 
acterized. Even though the measures are bipartite, it 
should be stressed that this is by no means an account 
of entanglement across different bipartite splittings. The 
maximally entangled state, with Ei = 1, Vi is the GHZ 
state [24 |. We have also analyzed some features of the 
classification induced by this set of measures and we have 
shown that the value of E\ determines whether a state 
is LU-equivalent to its complex conjugate or not. If this 
is the case, then four classes can be identified which are 
related to how the states may transform under determin- 
istic LOCC operations, showing further the physicality 



of our approach. It is worth remarking that such a con- 
nection can already be established from the value of a 
single measure. It will be interesting to study if the other 
measures can provide a finer classification of states un- 
der these transformations. Although the values of {Ei} 
depend on the chosen initial partition, it is worth re- 
marking that this classification is not (neither of course 
that of the states with the same entanglement). This 
reflects the fact that we are dealing with multipartite 
states and one might choose the partition depending on 
the particular task one wants to accomplish. Given the 
operational character of our approach, in the future we 
will study how to generalize it to more qubits and higher 
dimensions with the aim to understand the most rele- 
vant measures and identify possible new applications of 
multipartite entangled states. 

The research was funded by the Austrian Science Fund 
(FWF): Y535-N16 and F40-FoQus F4011-N16. 
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APPENDIX 

Appendix A: Mathematical properties of the 
3— qubit state decomposition 

In the main text we have proven that any 3-qubit state 
can be written up to LUs as 

|V) = -^=(|0>i|^>23 + \l)iU 2 ® Usl^hz), (5) 

where \ip s ) — a\00) + b\ll) with a > b has Schmidt de- 
composition and U 2 = Z(a)Y(j3)Z(~f) and U 3 = Y(/3') 
with Y, Z(£) — e^ (Ty -' respectively. Thus, any state is 
characterized by the five parameters, {E(tp s ), a, (3, 7, /?'} 
[H . In this section we prove that they can be taken such 
that E{ip s ) e [£(C 23 ),lE(C 2 a 3 )], a, 7 G (-7r/2,7r/2] and 
/3,/3' € [0, 7r/2] and show how the decomposition can be 
made unique. 

We start by proving that £(C 23 ) < E(ip g ) < ^(C^). 
As introduced in the main text it is useful to regard Eq. 
([5]) as providing a deterministic assistance protocol in 
which party 1 by performing some measurement will in- 
duce some entanglement between 2 and 3. More precisely, 
carrying out a Von Neumann measurement on the basis 
{|0), |1)} will surely produce 23 to hold a state with en- 
tanglement E(ip s ). The result then follows by taking into 
account that the concurrence and the concurrence of as- 
sistance, which are the extremal values (as measured by 
the concurrence) that can be obtained by any assistance 



protocol, can be obtained by some deterministic proto- 
col [HQ- It is not clear, however, that the particular 
protocol corresponding to decomposition ([5]) allows cer- 
tain states to attain any of these extremal values but this 
turns out to be the case (cf. Theorem 1 in Sec. IV below). 

To see that {a, /3, 7, /?'} can be taken to lie in the 
regions stated above, first notice that it suffices to 
consider the first and fourth quadrant since rotations 
around it correspond to the LU a\. The rest fol- 
lows by noticing that the following transformations on 
{a,/3,7,/3'} correspond to LUs: {a, — (3, 7, — (3'} (multi- 
plying by 1® cr z ®a z ), {a±7r/2, — /?, 7±7r/2, /?'} (because 
Y(p) = -Z{w/2)Y{-0)Z{v/2)) and {- 7 , — j9, -a, -/?'} 
(multiplying by a x ® Jj\ ® C/j). Notice that this im- 
plies that the transformation to {7, /?,«,/?'} (as well 
as {—a, /3, —7, /?'}) corresponds to complex conjugation, 
which is used in the main text. 

The question of uniqueness of the decomposition ([5]) is 
not crucial for our purposes as one can always define a 
particular choice of |-0 S ) , C/ 2 and U3 to make our measures 
{Ei} well defined. As we have shown in the main text this 
is enough to characterize operationally the classes of 3- 
qubit states with the same entanglement. Nevertheless, 
let us discuss this issue here for the sake of clarity and 
completeness. First, one may wonder if the choice of 
\ip s ) is unique (i. e. whether there is one single value 
such that E(\x Q )) — E(\xi))), This can be easily verified 
in almost all cases as one can read from Eq. (3) in the 
main text that oj is uniquely defined (up to the irrelevant 
addition of nir with n G Z) except when i) cq = Ci = 0, 
ii) pco = (1 — p)ci =^ and argc = 7r/2 and iii) c = in 
which it can take any value. Nevertheless, in case i) \ip s ) 
is also unique since C(|xo,i)) does not depend on u (in 
what follows in this case we will choose initially uj = so 
that the the unitaries U2 and U3 can be made unique). 
On the other hand, in cases ii) and iii) C(|a;o,i)) changes 
with oj and, hence, different choices for |?/> s ) are possible 
(which would lead to different values for our measures as 
Ei = E(ip s )). This is, for instance, the case of class 4 
inside the CLU class and, hence, it seems that there could 
be some physical meaning behind this fact. Nevertheless, 
the choice of \tp s ) can be fixed by taking the one for which 
E(ip s ) is maximal. Once the uniqueness of \ip s ) has been 
settled one could ask whether then /7 2 and U3 are unique 
(with the constraint that they must be respectively of 
the ZYZ and the Y form). In other words, whether there 
exist LUs {Vi} such that 

Vx ® V 2 ® U 3 -^(|0)i|V s ) 23 + \l)iU 2 ® U 3 \^ s ) 23 ) 
= -L(|0)x|^) 23 + \l)iU 2 ® U^ s ) 23 ). (6) 

It can be seen that this is possible iff Vi G {(Ji}t =0 up to 
a global phase Q. Hence, V 2 and V3 must be such that 
either V 2 ®V 3 \ip s ) 2? , = |^ s ) 23 and V 2 ®V 3 U 2 ® U 3 \i> s )23 = 
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U' 2 ® U^ia or V 2 ® 7 8 |V.)23 = ^2 ® ^1^)23 and F 2 ® 
IW2 ® ^3 1 ^3)23 = 1^3)23 (up to some possible change 
of signs). Then, if \ip s ) is not the maximally entangled 
state, since Z(£) (g> Z(—t;) for any £ is the most general 
LU operation that leaves |^ s ) invariant, it follows that 
the above considered transformations on {a, /?, 7, /3'} are 
the only ones that can occur. Hence, taking into account 
that the parameters are fixed to a, 7 € (— 7r/2,7r/2] and 
/8, j9' € [0, 7r/2] and imposing, for instance, \a\ > \j\ the 
decomposition is made unique. Notice, however that this 
is irrelevant as our measures {Ei} are invariant under 
these transformations. On the other hand, if \ip s ) is a 
maximally entangled state, then operations of the form 
U®U* for any unitary U are in this case the most general 
LU transformations that leave \if> B ) invariant. Then, in 
this case in order to make the decomposition unique and 
our measures {Ei} uniquely defined we take without loss 
of generality U 2 = Z{a)Y{(3) and U 3 = 1. 

Appendix B: Implementation of the controlled 
unitary operations to generate 5^ 

As explained in the main text, to apply the non-local 
gate C/ c 12 = |0)i(0|<K>l2 + |l)i(l|<8>C/2 on some state \4>) 12 , 
1 and 2 share in addition to their system qubits the CJ- 
state |*)i a i b2a2t = (|OO)i]0+> 2 + |H)it/ 2 ® l|0+) 2 )A/2 
3 (see Fig. 1 in the main text). Depending on the 
measurement outcome they obtain one of the four states 
Ul 2 (Ji ® <Jj\4>) (cto corresponding to the outcome <p + , a x 
to ip + , <Ty to ip~ and a z to 4>~). Thus, implement- 
ing U^ 2 and U^ 3 in this way leads to one of the 4 4 
states in the set {U^aj ® o 3 Ul 2 a\ ® erf |+)i |^ s ) 23 }. We 
show here that, due to the symmetry of the states, all 
of them are LU-equivalent to one of the four states in 
= {U^ 2 a 2 n \+) 1 \A) 23^=0- We have 

U 1 c 3 o-l®<T 3 U 1 c 2 v 1 k ®o- 2 \+) 1 \4> s ) 23 

~ LU U 1 c 3 U} 2 a 2 ®o- 3 m \+) 1 \^ s ) 23 
~ LU Ul 3 Ul 2 al\+) x \^ s ) 2 3, (7) 

for some m and n € {0, . . . , 3}. The reason for this is that 
commuting the Pauli operations with JJ\ % will at most 
change the sign of both, a and 7, and/or of (3 and/or 
of /?'. More precisely, the second line follows from the 
fact that a}\+) G {|+),|-}} and the fact that [7 c 13 crj 
is up to LU (from left) either U^ 3 or {U~l 3 )\ which is 
up to LU (from left) equivalent to U^ 3 a 3 . Note that 
the second case corresponds to changing the sign of f3' . 
The third line follows from the fact that the change of 
the sign of /3' can be accounted for by a change of the 
sign of p. Therefore, all output states are LU-equivalent 
to one of the states presented in Eq. (0. Now, since 
U™U™al\+)i\i) s ) 23 equals \ip) for n = 0, for n = 1, 
\ip*) for n — 2, and \ip') for n = z, this set coincides with 
the set S^. 



Note that contains at most four LU-inequivalent 
states, since under certain conditions some of these states 
are LU-equivalent (see main text). It can be easily ver- 
ified that the four possibilities are obtained with equal 
probability as for the implementation of U^ 2 

||ia c (0 + |2 b 2 c (0 + k!V^^a ; 2 =|*) lali) 2 a 2j+)lJV'o)2 c 3|| 2 

(8) 

always equals 1/16 and, hence, does not depend on 
i, j, fc, I and similarly for the posterior application of C/ c 13 . 

Appendix C: Identifying \ip) and \ip') with Eq 

We show here in detail that the binary measure Eq 
allows us to distinguish the two states \ip) and \ijj'}, 
if they are not LU-equivalent. Recall that E(j(\&)) 
is defined as follows: E 6 {\^)) = if E 1{23 {\^}) = 
min|$ i ) e 5^J5 1 |23(|$i)) and E 6 (\$)) = 1 otherwise. Since 
£i|23(l$)) = £i|23(l$*)), for any state |$), the set 
{Ei\ 23 (\&i)), l$i) € S^} contains at most two values. 
Using the same reasoning as in the main text, £d|23 is 
uniquely given by the value of 

\(i/>.\Z(a)Y(0)Z(r,) ® Y{p')\^ s )\ 2 = 4a 2 & 2 sin 2 /W / 3' 
cos 2 (a - 7) + cos 2 cos 2 P'[{a 2 - b 2 ) 2 + 4a 2 6 2 cos 2 (a + 7)] 
+ 4a6 cos P cos p' sin/3sin/3' cos(a + 7) cos(a — 7). (9) 

Recall that if \ip) corresponds to the parame- 
ters (E(\ip s )), a, P, 7, P') (see Eq. (2) in the main 
text), then \ip') corresponds to the parameters 
(E(\ijj s )), a, — P, 7, P'). Thus, the first two terms in Eq. 
(J9j) take the same value for both \tp) and l^'), while 
the last one differs in sign. Hence, E^ 23 takes dif- 
ferent values for these states and therefore Eq distin- 
guishes them, unless the last term in Eq. Q is zero. 
We show in the following that this is the case iff \ip) 
and \ifj') are LU-equivalent. If p, p' = 0,7r/2 \ip) and 
\ip') are clearly LU-equivalent. If a + 7 = tc/2 or 
a — 7 = 7r/2 they are LU-equivalent, as can be seen by 
using that Y(p)Z(ir/2) = Z(n/2)Y(-0) together with 
the LU-equivalences of {a, P, 7, P'} given in Sec. I. It re- 
mains to consider the case a = or b = 0. Then, \tp) 
is LU-equivalent to ^/p|0> |00> + y/T^p\l)U 2 (a, P, 7) <g> 
Y(P')\00), which is obviously LU-equivalent to \ip'), since 
a z \0) = |0) and Y{p')a z = a z Y{~P'). 

Appendix D: Characterization of CLU and NCLU 
classes 

The aim of this section is to prove the characterization 
of the CLU class. We will first derive a mathematical 
characterization (see Lemma I below), which will then 
be used to classify CLU states according to the value the 
measure E\ takes (Theorem 1 below). 
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As in the main text we consider without loss of gener- 
ality the following decomposition of 

iv) = v^|o}iVo> + \/r r P|i>i^i), (io) 

where \ipo) — ax|00) +a4|ll), with a,i,a± £ M. and |V>i) = 
6i|00) + 6 2 |01) + 63 1 10) + 6 4 |11), with b, G C. Moreover 
the following conditions are satisfied: 

ci) (VilVi) = kr 

C2) co,d > 0. 

Note that such a decomposition is always possible since 
a local phase gate, diag(e lQ °/ 2 , e Ml / 2 ) on the first qubit 
changes Co (ci) to c e M ° (c\e lai ) respectively. Given this 
decomposition c = —a\b^ — a^bi- We are going to show 
now that c is real or purely imaginary iff |-0) is CLU. 

Lemma 1 For a state \ip) there exists a decomposition 
as in Eq. U7jj) with c £ R or c £ iR iff is CLU. 

Proof. First, notice that it suffices to prove that the state 
is CLU iff there exists a decomposition with cq,ci,c £ 
R. This is because the latter condition holds iff there 
exists a decomposition with cq, c\ > and c £ R, iR as a 
local phase transformation Ck Cke lT! (k = or k = 1 ) 
induces c — > ce !7r / 2 . Also, we note that the condition that 

= (ipoltpi) = aib% + 0464 implies that the phase of 61 
equals the one of 64 (unless cq = 0, which will be studied 
separately). We write b\ = b\e 1 ^ and 64 = b\e l<t) ', where 
both, b\ and b\ are real. Let us start by proving the "only 
if part. That is, we assume that c = —0164 — 0461 is real 
and show that this implies that \ip) is CLU. c being real 
implies that either (j) = kit for some integer k or c = 0. 
In the first case, c\ = 2(—b r 1 b r A + & 2 &3), which is real iff 
62^3 is real. Hence, b 2 = b\e %x and 63 = b^e~ lx , where 
b\ 3 are real. It is then easy to see that the local unitary 

1 ® Z(—\/2) <£> Z(x/2) transforms into a state with 
real coefficients, which proves that the state is CLU. If 
c = — ai&4 — a^bi = and a\ ^ a\, CI further implies 
that 61 = 64 = . Then, it can be easily shown again 
that the state can be mapped by local phase gates to a 
state with real coefficients, which is CLU. If a 2 = a 2 , then 
without loss of generality we can take a\ = a 4 = l/v2) 
i. e. |-0o) = \<t> + )- In this case CI does not add any 
further constraint and we have 64 = —61, i. e. = 
61 1 00) + 6 2 |01) + 63 1 10) - &i|ll>- We will show again 
that there exists local unitaries that transform \ip) into a 
state with real coefficients. Let X = &i|0)(0| + fo 2 1 0) < 1 1 + 
63 1 1> <0| - &i|l)(l| and let X = UEV^ be its singular 
value decomposition. This means that the local unitaries 
W and V T transform into its Schmidt form, i. e. 
W (g> U T |?Ai) = IV'f) £ R 4 . Notice that the fact that 
TtX = implies that Tr(T,V^U) = and, hence, V^U = 
e^W, where W is a unitary matrix with real diagonal 
entries i»u and W22 and off-diagonal entries of the form 
w%2 = w 12 e 1 ^ and W21 = u> 21 e - ^ with w\ 2 and w 21 both 



real. Let |$) = 1 ® Z{-a)U* ® Z{a)V T \ip). Using that 
the local unitary transformation Z{a) ® Z(—a) for any a 
leaves any state in the Schmidt form invariant and that 
U\ ® U 2 \(j) + ) = 1® U 2 Ui\<p + ) for all unitaries [/ 1>2) we 
have that 

= ^e-^\Q)l®Z{a)W*Z{-a)\<j ) + ) + ^r^p\l)\i,{). 

(11) 

Taking into account the above conditions on W, we have 
that Z{a)W* Z(—a) is a real matrix if a is chosen such 
that a — £/2. Therefore, the local unitary diag(e' 1 '^, 1) <E) 
Z(-£/2)W ® Z(£/2)V T transforms into a state with 
real coefficients as we wanted to show. 

It remains to consider the case cq — 0. We choose 
without loss of generality 04 = 0. Then, condition (CI) 
implies that b\ = and therefore c = —64, which is real 
iff 64 is. Condition C2 implies that 6 2 &3 is real. In this 
case the local unitary 1 ® Z(— x/2) ® Z(\/2) for some 
proper choice of x transforms |t/>) into a state with real 
coefficients, which completes the proof of the "only if 
part. 

Let us now show the implication in the opposite direc- 
tion. If |V>) is CLU, it has been shown in [5[ that, then, 
there exists a product basis in which the state has real 
coefficients, i. e. \ip) ~ L u J2ijk ^ijk\ijk) with X ijk £ R 
V i,j,k. Thus, the singular value decomposition of the 
(real) matrix X — J2ijk ^ijk\i){jk\ is real and, therefore, 
so is the Schmidt decomposition of the state in the split- 
ting 1 1 23. Hence, cq, c\ and c can be chosen to be real. 

□ 

We are now in the position to prove that the CLU 
class can be characterized via the values the entangle- 
ment measure E\ = E(\ip s )) takes. Let us recall here the 
decomposition [Eq. (2) in the main text] we have proven 
in the main text, as we will be using it throughout this 
Supplemental Material. Any 3-qubit state can be written 
up to LUs as 

|V) = -^(|0) 1 |V.) a8 + \lhU2 ® CTahMaa), ( 12 ) 

where \ip a ) has Schmidt decomposition, [/ 2 = 
Z(a)Y(P)Z(r,), U 3 = Y(j3>), Ety.) £ [g(C 23 ), £(C 2 a 3 )], 
a, 7 G (-7r/2,7r/2] and 0,0' £ [0,tt/2]. 

Theorem 1 £ CLU iff there exists a decomposition 
of the form such that either E\ = (£(C 2 3) or E\ = 

£(c 2 a 3 )- 

Proof. First, it should be noticed that, since C{p2z) = 
si - s 2 d and C a { P 2z) = s\ + s 2 |4j] where {st} are the 
singular values of r (arranged as usual in non-increasing 
order), it can be easily seen that C 2 3 = C_ and C^ = C + 
with 

C 2 ± =p 2 cl + {l-p) 2 cl + 2p(\-p){\~c\ 2 ±\c QCl -c 2 \). (13) 
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On the other hand, we have that 

C(^) 2 = p 2 cl + (1 -p) 2 c 2 + 2p(l -p)(coCi cos 2w + 2|6| 2 ), 

(14) 

with lu given by Eq. (3) in the main text. From this, it is 
then clear that if c £ M. or c £ iM., then either C(V's) = C+ 
or C(ips) — C-, which together with Lemma 1 proves the 
only i/part of the theorem. To see that the implication in 
the opposite direction holds as well, we have to prove that 
argc = mr, nn/2 are the only solutions to C(tp s ) = C± 
under the constraint (3) of the main text. This constraint 
can be rewritten as 



cos 2lo 



1 



(1 



i cos(2 arg c) 



1 



1 



t cos (2 arg c) 



(15) 



where we use the notation x = (pco + (l — p)c\)/{pcQ — (1 — 
p)c\). Using this and writing y = cos(2argc), C{i)j s ) = 
C± boils down to 



C Ci 



x 



(l + x 2 )y + 



1 + x 2 — (1 — x 2 )y 
= ± 



44 + \c\- 



2c ci\c\ 2 y. 



(16) 



After squaring this last expression, our problem reduces 
to identifying the zeros of a third order polynomial p{y) ■ 
More precisely, we know two zeros y\ = — 1 and 2/2 = 1 
(which correspond to argc = mr 1 mr/2) and we have to 
check that the third one 2/3 fulfills 2/3 ^ (—1, 1), so that 
t/3 = cos(2 argc) cannot be inverted and, hence, there are 
not more solutions to C(ip s ) = C±. Indeed, using that 
|p(0)| = I72/12/22/3I where 7 is the leading coefficient of p, 
we see that 



m\ 



|2c oCl x 2 + |c| 2 (x 4 -l)| 
|c| 2 (x 2 - l) 2 



(17) 



Now, since we take co,ci > (which also implies that 
x > 1), it is easily seen that \y^\ > 1, which completes 
the proof. 

□ 



Appendix E: Characterization of the 4 subclasses 
inside the CLU class 

In this section we will prove the characterization of the 
four subclasses in which the CLU class can be divided ac- 
cording to the particular value of E\ (which by Theorem 
1 above can be either (£(6*23) or €(Cf 3 )). 

The fact that class 1 (Ex = £(C 23 ) = £(C 2 Q 3 )) corre- 
sponds to the W class can be easily seen as follows. It can 
be shown that (C^) 2 = C| 3 + T123 (see e. g. Q), where 
T123 is the tangle. Since the W class is characterized by 
a vanishing tangle, the result follows straightforwardly. 

The identification of the other three classes will re- 
quire more work. As Lemma 2 below shows, the sign of 



det r will play a crucial role in this task. We recall that 
the matrix r associated with a state \ip) is a function of 
its reduced state p 2 3- It is given by Ty = JPiPji^i^j) , 
where {p{\ are the eigenvalues of P23 and {|"0i)} the corre- 
sponding eigenstates. Therefore, it is symmetric. Since 
detr = iff the tangle T123 = (as follows from 0]), 
which corresponds to class 1, we assume in the following 
that detr 7^ 0. First note, that by redefining the sign 
of Ci £ R we can always obtain a decomposition with 
c G K (recall that we are dealing with the CLU class 
and, hence, Lemma 1 applies). Thus, r £ R 2x2 in what 
follows unless otherwise stated. 

Before relating the sign of det r to our classes, we need 
to discuss when this property is well defined. This is 
because it could in principle happen that two LU equiv- 
alent states lead to different signs for detr. We now 
prove that this is indeed the case iff there is a choice of 
eigenstates \tpo) and \ipi) such that c = 0. First, if the 
latter condition holds, then it is clear that under a dif- 
ferent choice of global phase for \ipo), \ipo) — !• i\ipo), r 
remains real and the sign of its determinant is flipped. 
To prove the implication in the other direction we have 
to consider how r transforms under the LU transfor- 
mations \ip) -> W\ W2 ® W 3 |V')- Using that t = 
M^a y ®a y M\, where My = ^/p\0) (Vol + vT-~p|l) (M 
and that Wa y W T = det(W)a y for any matrix W, we 
find that t transforms to 



det(W 2 ) det(W 3 )WiTiy 1 T . 



(18) 



Now, we have to distinguish the following two cases, 
i) p / 1/2 and, hence, the eigenspace of P23 is non- 
degenerate and the Schmidt form unique. Then, W\ 
can only take the form W\ = diag(e lai ,e la2 ) to pre- 
serve the orthogonality of \ip ) and |^i). In this case 
det f = - det r iff det(I^ 2 ) det(VU 3 )e l(Ql+a2) = ±i. How- 
ever, notice that for the off-diagonal entry we have that 
t 12 = det(iy 2 )det(WAs)e i(Ql+Q2) Ti2 and, hence, the fact 
that both t and f should be real implies that c = has 
to hold, ii) p = 1/2. Then, the eigenspace of P23 is de- 
generate and the Schmidt form not unique (as \ipo) an d 
\ipi) can be chosen to be any orthonormal basis of the 
eigenspace). Thus, W\ might be an arbitrary unitary 
matrix. We show now that in this case there always ex- 
ist LUs W\ <g> 1^2 <8> W3 such that the sign of det r can 
be changed and such that c = 0. Let V denote the real 
orthogonal matrix that diagonalizes r and let us choose 

Wi = ( * ° J V and W 2 = W 3 = 1 in Eq. (JTSJ). The 

matrix W\ defined in this way is such that det W\ = ±i 
(thus det f = — det r) and with this choice f is real and, 
furthermore, diagonal. Hence, there exists a choice for 
which 2 = 0. This ends the proof. 

Therefore, special care must be taken with those states 
that allow for c = 0. This is not surprising, as we have 
seen in Sec. I that these are the only ones for which the 
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choice oi\ip s ) in Eq. (|12j) is in principle not unique. Let us 
now show how the sign of det r is related to the different 
classes. 

Lemma 2 There exists a decomposition of the form 
such that Ei — <£(C23) iff det r > and such that E\ — 
<£(C$ 3 ) iff detr <0. 

Proof. Since r g R 2x2 , the unitary in Eq. (1) in the main 
text that allows for the transformation to the decompo- 
sition (TT2"]) can be taken to be U = y j /\/2. We 

will also use that | det r| = T123/4, which follows from the 
results of Q. With this in mind, 

W4= I dot r I - \dct(UrU T )\ = i|C? - (Trr) 2 |, (19) 

where C\ is such that E\ ~ l£(Ci). On the other hand, 
since r is symmetric (and now real), it is Hcrmitian and, 
hence, the absolute value of its eigenvalues equal the sin- 
gular values. Therefore, either |Trr| = s\ — s 2 — C'23 
(when detr < 0) or |Trr| = si + s 2 = C 23 (when 



detr > 0). Now, since (C| 



a 2 

°23 



T123 has to hold 



in Eq. (fll?|) , the result follows. (This can also be seen by 
comparing Eqs. HU) and (fTI])). □ 
To identify the sign of det r we will make intensive 
use of the standard form (up to LUs) for 3-qubit states 
presented in 

\4>) = Ao|000) + A ie ^|100) + A 2 |101) + A 3 |110) + A 4 |lll), 

(20) 

where A^ > Vi and <fi £ [0, w]. Some equations might 
be more compactly written in terms of the polynomial 
invariants [j| 

Ji = |A 1 A 4 e < *- A 2 A 3 | 2 , 
J 2 = AqA 2 , 
•h = AqA 2 , 

J4 = AqA|, 

Js = A (Ji + A 2 A 2 - A 2 A 2 ). 



(21) 



Grassl's invariant, which discriminates \ip) and \ip*) for 
states in the NCLU class, is given by 

J 6 = XtXl(\4(l - 2A 2 - 2Xl) + 2A 1 A 2 A 3 e^) 2 . (22) 

According to @, the CLU class can be characterized by 
certain conditions satisfied by the polynomial invariants 
[Eq. (|2ip]. Namely, a state is in the CLU class iff either 



\M = ^Ji^Jz 



(23) 



(J 4 + J 5 ) 2 - 4(Ji + J 4 )(J 2 + + Ji) = 0. (24) 



The states satisfying Eq. ([23)) are moreover shown to be 
equivalent to those for which the standard form given in 



Eq. ([20]) is already real (i. e. <p = 0,tt) [5]. Alternatively, 
these two subclasses can be characterized by the value 
of Jq. Since they are CLU, Je{ip) = J&{^*) and it fol- 
lows that Jq must be real. The subclass characterized by 
Eq. (|23p corresponds to Jg > (the expression which is 
squared in Eq. ([2"2"j) is real) and the subclass characterized 
by Eq. ([24]) to Jq < (the expression which is squared in 
Eq. (|22p has to be purely imaginary). Furthermore, it is 
shown in Q that the states for which Eq. ([23]) holds can 
be written up to LUs as 



ci|000) +c 2 \(j) 



19293/ 



(25) 



with 



Ci G 



and 



which correspond to our class 



2. On the other hand it is also shown therein that the 
states satisfying Eq. (pM)) are LU equivalent to states of 
the form 



|OOO)+e 45 |<^20 3 ) 



(26) 



up to normalization, where S £ R and ^, e I 2 . Notice 
that this class corresponds to our class 3 (loj . 

Our aim is to connect these families with the sign of 
det t. More precisely, we will show in the following that 
for states for which Eq. (|23| holds, i. e. they can be writ- 
ten as in Eq. (1231) . (Eq. (IM1) holds, i. e. they can be 
written as in Eq. ^) then detr < (detr > 0). Us- 
ing then Lemma 2, we have that any state in class 2 
(class 3) fulfills E x = (E(C| 3 ) (E 1 = £(C 23 )). Since con- 
ditions (|23[) -(j24 j) characterize completely the CLU class, 
the implication has to hold as well in the opposite direc- 
tion, i. e. if E x = £(Cf 3 ) (£1 = €(C 23 )) then the state 
can be written as in Eq. fl5SJ| (Eq. (|26]l). Thus, this will 
prove the characterization of classes 2 and 3. Of course, 
this makes sense only if the sign of det r is well defined. 
Hence, for this task we will assume that a choice of r for 
which 5 = does not exist. We will see afterwards that 
this assumption is justified as we will prove that c can be 
taken to be zero iff both conditions (|23[) and ([24]) hold at 
the same time. This will finally lead to the characteriza- 
tion of class 4 which is given by the intersection of classes 
2 and 3. 

First, we need to write r in terms of the A's. Hence, 
we need to find the 1 1 23 Schmidt decomposition of the 
state given in Eq. ([2"U]) . Notice that Lemma 2 holds when 
the global phases of the Schmidt eigenstates for 23, |i/> ) 
and are chosen such that r £ R 2x2 (which is always 
possible for CLU states by Lemma 1). Thus, once r is 
computed from the Schmidt decomposition of the state 
of Eq. ([20]) a further transformation might be needed for 
the latter condition to hold. As before, we have that 
these global phase transformations transform the entries 
of t as 



c -> c e 



Ci —> C\e 



iff 



de i(a+f3)/2_ 



(27) 



We will denote by t t the real matrix obtained by these 
transformations. However, fortunately, since, as men- 
tioned above, for the subclass characterized by Eq. ([2"3"]) 
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it holds that (f> = 0, 7T, we will always deal with real ex- 
pressions in this case. Hence, we will only need to take 
care of this fact for the second subclass of states satisfy- 
ing Eq. 423}. 

Let us start by considering the case in which Eq. ([23]) 
holds and Eq. ([24]) does not. We have to show that this 
implies that detr r < 0. First, we consider that Ai = 
(which agrees with the fact that Eq. ([2"3"]) is satisfied), 
as in this case the state in Eq. ([20"]) has already the 
Schmidt form. Thus, p = A 2 ,, cq = 0, c\ = 2A2A 3 /(1 — p) 
and c = — A 4 /\/l — p. Thus, r is then already real and 
det t> = — J 4 < as we wanted to show. Notice that here 
one should exclude A 2 , = 1/2, as this would imply using 
Eqs. ([ST) that Eq. ([24]) would hold as well. This is not 
surprising since this corresponds to p = 1/2 for which we 
have seen that the sign of det r r can be changed by LU 
transformations. We will see shortly that the Ai = case 
is the only possibility in which p = 1/2 may hold. 

We now consider the case when Ai 7^ 0. To compute 
the Schmidt form, we need to find the singular value 
decomposition X = ITEV^ of 



and, moreover, 



X 



A 

Aie ,; ^ A 2 A 3 A 4 



(28) 



U providing the Schmidt basis for qubit 1, the singular 
values <t± providing the Schmidt coefficients (i. e. 0+ = 
p, <r_ = 1 — p) and the first two columns of V providing 
the Schmidt basis for qubits 2 and 3 (i. e. v\ — |V>o) and 
t'2 = I ipi ) ) ■ After some algebra one finds that 



ct± 



l±y/l- 4( J 2 + J 3 + J 4 ) 



U = 



(X 2 -a+)/k, 



-AoAie^/fc- 
(A 2 - <r_)/fc_ 



(29) 



where k± are just normalization factors to make the 
columns of U of unit norm (i. e. k± = A 2 A 2 + (A 2 , — cr±) 2 )- 
Notice that the case a± = 1/2 does not occur here 
as J2 + J3 + J4 = 1/4 implies, using Eqs. ([STjl . that 
Aq(1 - A? - A 2 .) = 1/4, and therefore that Ai = 0. As 
VYJ = X^U, we then have that 



/ -A^Aie^/fci + (A 2 - a ± )X ie ^/k ± \ 
A 2 (A 2 - <r±)/fc± 
A 3 (A§ - a ± )/k ± 
\ 4 (\q - cr±)/k± 



a±v± 



(30) 



where v\ — 14. and i> 2 



From this, it follows that 



fc±(T±co,i =2AgAiA 4 (Ag-(7 ± )e i * 
-2(X 2 -a ± ) 2 (X 1 X 4 e^- 
k+k^^/d+FZc = A^AiA 4 (2Ag - l)e^ 

+ 2AgA?(AiA 4 e < *-A 2 A 3 ) 



A2A3), 



(31) 



klk 2 



2 ;„2 _ _ -2 



k\kia + a-c z = A^A^A^(4A^ - 4A ( 
where 



1 e 



2i<p 



A, (32) 



a = 4x^x1(2x1 + xj-r 



+ 4A4a 1 A 2 A 3 [A 4 (1 - 2X1 ~ 2A?)e** + AiA 2 A 3 ]. (33) 



One then finally finds that 



klk 2 _ det r = AX^XlXji J 2 + J 3 + J 4 



(34) 



As discussed above, then for the states for which Eq. 
(|23| is true it holds that k 2 + k 2 _ det r r = 4AgAf A 4 ( J 2 + J3 + 
Ja — 1/4). Since it can be shown that J 2 + J3 + J 4 < 1/4 
[5| , we have that det r r < in this case. This finishes the 
identification of class 2. 

Let us now consider the case for which Eq. ([2"3"]) does 
not hold (i. e. <j> ^ 0, rr) but for which Eq. ([24]) docs 
Hence, we need to prove that detr,- > 0. Notice 
that Ai = implies that Eq. ([23]) is true. Therefore, we 
just need to consider the Ai 7^ case which, as above, 
leads us to Eq. (jM]) . However, as discussed before, to 
construct r r for this subclass, a specific transformation 
of the form l|27p must be performed. That is, we choose 
the global phase of |^o) [ca/2) and the global phase of 
I -01 ) (/3/2) such that r gets real. Since under such a 
transformation detr changes to detTe^ Q+/3 \ we get a 
real expression in ([34]) iff either a + j3 = —2(f) (which 
would lead to det r < 0) or a + f3 — ir — 2<fi (which would 
lead to detr > 0). Moreover, ce*( a+ ^/ 2 must be real, 
i.e. c 2 e l ( a+ ^ > must hold. In other words, the sign of 
g2 e -2i0 tell us whether a and {3 have to be chosen 
such that a + /3 = —2<f> or a + (3 — ir — 2(j). In the following 
we prove that for any choice of the phases a, /3 which lead 
tore K 2x2 , it holds that c 2 e~ 2 ^ < 0. This means that 
a+f3 = ir— 2(f> is the proper choice to obtain r r . Hence, for 
the CLU states such that (f> ^ 0, ir (i. e. those satisfying 
Eq. ([2"4"]0 . it holds that detry > and therefore these 
states belong to class 3. To see that c 2 e~ 2 ^ < 0, notice 
that the transformation ([2"T]) must fulfill both cqCi E M 
and c 2 el [IH ■ Equation (|32l) implies that for the states 
for which this is possible Ae~ 2% ^ £ R must hold. This 
means that these states are such that 



A 4 (l - 2Aq - 2A 2 ) sin <j> + X 1 X 2 X 3 sin 20 = 0. 
Since now (f> 7^ 0, rr, Eq. (1351) is equivalent to 

A 4 (l - 2A 2 , - 2A 2 ) + 2AiA 2 A 3 cos0 = 
or, simply 

A 4 (l - 2A 2 - 2Xj) 



COS( 



2AiA 2 A3 



(35) 



(36) 



(37) 
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Inserting Eq. (J55J in Eq. ([32J leads to 

k 2 + k 2 _ a + a-d 2 e- 21 * = 4A4A 2 A 2 [l/4 - {X 2 + A 2 )(l - X 2 - A?) 

+ 4A^A?A 2 A 3 [A 4 (1 - 2Xl - 2A 2 ) cos0 + AiA 2 A 3 cos 20]. 

(38) 

The second line of the above equation can be rewritten 
using Eq. (|36|). which leads to 

A 4 (l-2A 2 l -2A 2 )cos0 + AiA 2 A 3 cos20 = -AiA 2 A 3 . (39) 

Using now that 1/4- (Ag + Af)(l - - A?) < (1 - 2A§ - 
2A 2 ) 2 for arbitrary Ao, Ai such that A 2 , + X\ < 1 we find 
the following upper bound for the left hand side of Eq. 
©8) 

< A^A 2 [A 2 (1 - 2A 2 - 2A 2 ) 2 - 4A 2 A 2 A 2 ] < 0, (40) 

where in the second inequality we have used Eq. (1371) . 
Thus, c 2 e~ 2l< ^ < 0, which concludes the proof that for 
any state which does not fulfill Eq. (|23| but does fulfill 
Eq. ([24]) det t > holds. 

Finally, let us consider the case for which both Eqs. 
(f2"3")) and ([H]) hold, i. e. J 6 = 0. Recalling that we are 
excluding the class of W states already characterized as 
class 1 (for which these conditions are also true), this 
means that the states are LU equivalent to 

|000) + |0i020 3 ) (41) 

up to normalization and with 0j £l 2 . This can be un- 
derstood from the fact that no state in class 2 can be 
transformed to any in class 3 via LOCC. In particular, it 
cannot be LU equivalent to any state in class 3. There- 
fore, this class must contain the intersection of the forms 
given in Eqs. (I25[) and (|26p . which corresponds to our 
class 4. It remains to show that this class is the only 
one for which detiv can take both signs. That is, leaving 
aside the W class, class 4 is characterized by the pos- 
sibility of having c = 0. We will prove that both Eqs. 
(f2"3")) and (f2~4"|) hold iff there exists a choice of r r for which 
c = 0. The proof of the if part is immediate as can be 
seen as follows. If p = 1/2 we have seen before that this 
already implies that c = is possible. If p ^ 1/2 ,we have 
just seen that if Eq. (f23|) holds (0 = 0, n) the matrix r 
constructed from the standard form ([20]) is already real 
and that c 2 > 0. If Eq. (JM]) is moreover satisfied, then 
it also holds that c 2 < and, hence, c = 0. To prove the 
implication in the other direction we assume that 2=0. 
If Ai 0, we can use Eq. (|3"Tj) . which then yields that 

A 2 [A 4 (1 - 2A 2 - 2A 2 ) + 2AiA 2 A 3 e-^] = 0. (42) 



This means that J 6 = and, thus, both Eqs. (|23|) - (|241) 
hold. Following our discussion before Eq. (|28|) . in the 
] case Ai = 0, there exists a choice of r r with c = iff 
either one of the two possibilities is fulfilled: A 4 = or 
p = 1/2, for which the Schmidt form is not unique. The 
first one clearly leads to Jq = (and, in particular, to 
the W class since J 4 = Ti 23 /4 [5|). In the second case, 
Aq = 1/2 must hold as well and, therefore, J§ = 0. This 
completes the proof. 



Thus, as we discussed in Sec. I, since c = 0, for the 
states in class 4 there exist in principle different choices 
for \ijj s ) in the decomposition (TT2J) . 
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